Kinetics of self-induced aggregation in Brownian particles. 
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We study a model of interacting random walkers that proposes a simple mechanism for the 
emergence of cooperation in group of individuals. Each individual, represented by a Brownian 
particle, experiences an interaction produced by the local unbalance in the spatial distribution of 
the other individuals. This interaction results in a nonlinear velocity driving the particle trajectories 
in the direction of the nearest more crowded regions; the competition among different aggregating 
centers generates nontrivial dynamical regimes. Our simulations show that for sufficiently low 
randomness, the system evolves through a coalescence behavior characterized by clusters of particles 
growing with a power law in time. In addition, the typical scaling properties of the general theory 
of stochastic aggregation processes are verified. 
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I. INTRODUCTION 

A process frequently encountered in the study of chem- 
ical and physical phenomena is the aggregation of small 
particles joining each other to form larger spatial struc- 
tures and clusters. The comprehension of the gen- 
eral properties of aggregation processes constitutes a 
cross-disciplinary interest for pure and applied research 
[E 0, S 0, IE @| with broad implications to engineering 
and industrial technology Q. 

Aggregation is also a basic process in those biologi- 
cal systems where cooperation activity among individ- 
uals usually involves social behaviors. A well known 
example is animal grouping where an ensemble of indi- 
viduals belonging to the same species live together into 
organized communities such as insect swarms, mammal 
herds, fish-schools, and bird flocks [1, Q. The emer- 
gence of cooperation is still one of the most puzzling 
mechanism occurring in biology, mainly because social 
and altruistic behaviors are against evolutionary selec- 
tion which, instead, promotes antagonism and competi- 
tion. Mathematical population biology [IE HH, [HJ and 
Game Theory [IE HE HE HE] have always attempted 
to reconcile altruism with evolutionary selection trying 
to establish rigorous and quantitative basis to the emer- 
gence of cooperation among individuals. One of the first 
approaches, known as kinship theory, has been elaborated 
by W.D. Hamilton [l7|. It is based on the principle of 
kin selection and allows altruism to arise among siblings 
provided they share enough common genes. Although 
the Hamilton's principle strictly works for related indi- 
viduals, it has been also applied to cases of absence of kin 
recognition as a rather general explanation for altruism. 
Alternatively, when genetic arguments cannot be mani- 
festly invoked because cooperation involves unrelated in- 
dividuals, it is reasonable to resort to minimalist phe- 
nomenological models describing the mutual advantage 



to form a group or to make coordinated movements [181 ] 
in terms of effective "social" interactions. Synchronized 
and collective behaviors, on the other hand, are the result 
of "short-" and "long-range" interactions between differ- 
ent units which, in social biological systems, are not very 
dissimilar to those used in the description of fluids and 
condensed matter. In this perspective, it is reasonable to 
hope that some general features of collective behaviors of 
living organisms can be simulated and classified through 
the same principles of statistical mechanics successfully 
applied to understand the structural organization of the 
matter. 

General and simple theoretical descriptions of cooper- 
ation in social communities requires two principal ingre- 
dients: i) the space of states, where each point represents 
the status of the individuals and ii) the strategy, that is 
the rule according to which the individuals (players) de- 
cide to change their status in response to partial or com- 
plete information about the actions of the other players 
[161 [l^ | . The prevalence of cooperation is the outcome 
of those strategies rewarding altruistic acts and punish- 
ing defections. These rules should be given as simple 
and generic as possible in order to capture the essential 
and universal aspects of the problem without making the 
theoretical approach extremely complex or too computa- 
tionally demanding. 

In a previous work [2(| HH we studied a model sug- 
gested by Sigmund and Nowak [22j to test the hypothe- 
sis of the emergence of cooperation by indirect reciprocity 
among unrelated individuals. In their game theoretical 
approach, Sigmund and Nowak assumed that cooperation 
to work in evolved social systems requires the knowl- 
edge of "reputation or status" of their members (play- 
ers). Thus a dynamical coordinate, S, the image score, 
is assigned to each player signaling his/her reputation or 
status to the group. S is updated according to the al- 
truistic or selfish acts made by players in the past and it 
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is an information accessible to the whole community. At 
each turn of the game, a randomly selected player with 
score S assesses the possibility to provide help only to 
those opponents with a score greater than S. The image 
score, thus, determines the selfish or altruistic strategy. 
Altruism, despite the cost, is preferably since it increases 
S that, conversely, is lowered by a defection. In this 
scheme, cooperation occurs thanks to the principle "help 
and you will be helped" . 

A simpler version of this problem was recast, in a re- 
cent paper [20], in terms of a non- linear Fokker-Planck 
equation for the population of individuals with a certain 
image score. The equation had a non-local drift term 
that characterized the strategy. The peers exert a sort of 
"pressure" on each other within a finite range of influence 
in order to uniform their image score to the majority. In 
this sense, the model produces an aggregation mecha- 
nism in the space of scores driven by a population gra- 
dient; aggregated states correspond to situations where 
cooperation is achieved. The model, upon changing the 
system parameters, exhibited a transition from an aggre- 
gation behavior to a uniform state with no prevalence of 
selected image scoring. 

This model, in another context and with proper 
changes, can also be used to describe pattern forma- 
tion and chemotaxis phenomena where diffusion com- 
petes with a drift induced by chemical or population gra- 
dients H,[H,[2E[2E[i3. 

In this paper we focus on the dynamics of th e ag gre- 
gation process which was not considered in Ref . [20| • To 
this aim, we employ a discrete, or individual based, vari- 
ant of that model, where a set of individuals modify their 
status according to stochastic differential equations cor- 
responding to a Brownian motion with a drift induced by 
the spatial distribution of other walkers. This model can 
also be interpreted as a collection of interacting random 
walks [28] where the path coalescence depends on the bal- 
ance between the diffusion and the nonlocal drift. The 
particle-based modeling permits one to track the state 
and position of each component of the system. It is ap- 
propriate to study the fluctuations in this kind of ag- 
gregation dynamics characterized by the formation and 
merging of growing clusters. We shall see that, differently 
from previous proposed mechanisms, the coalescence of 
trajectories does not emerge as a result of direct inter- 
actions between particles [29j or via reaction-diffusion 
mechanism [8j, but rather it is the consequence of the 
drift term sensible to density fluctuations occurring even 
at distances relatively large. We found that in this model, 
aggregation and clustering follow the general sca ling be- 
havior of stochastic coalescence phenomena [3(| HJ [32| 
with power law for the kinetics of the number and average 
mass of clusters. 

The paper is organized as follows. In sect. II we de- 
scribe the model and the basic features of its dynamics. 
Sect. Ill is devoted to the presentation and discussion of 
numerical results. Finally, in sect. IV, we draw conclu- 
sions. 



II. MODEL OF INTERACTING INDIVIDUALS 

The model consists in a system of N units (individuals) 
which change their state Xi according to a sort of major- 
ity rule. The variable Xi, referred to the i-th member, 
might indicate the reputation score in indirect cooper- 
ation models, the position in a possible chemotaxis de- 
scription, or some other amplitude characterizing the role 
of an individual within a population biology framework. 
Hereafter, without loss of generality, it is convenient to 
adopt the terminology of spatially distributed systems, 
thus x will be a spatial coordinate. The formulation, 
however, can apply to other and different contexts. 

We assume that each individual changes Xi according 
to the non-local stochastic equation of motion 



Xi = v(xi) + V2D& 
with the drift v defined by the formula 

w + (xi,t) - w-(x h t) 



v(x l ) = A 



w + (x,i,t) +w-(xi,t) 



(1) 



(2) 



where uu± are given by 



w±(xi) = ^2e[±{xj 
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Xi)] exp(-a\xj - xt\) (3) 



The velocity at which an individual "decides" to move to 
the left or to the right is the result of his perception of the 
unbalance, u;_|_ — w_, between the populations at his left 
and right. This perception is simply modeled by the ex- 
ponential weight with a coefficient a = 1 /rg that defines 
the sensing distance, i.e. the range ro within which one 
individual still perceives the presence and the influence 



of the other members of the group 



<r ). This 



means that, unless a = 0, only partial information about 
the population is available to each individual. In this 
scheme, Vi is likely to point towards the most populated 
regions if not too far from the i-th player; the evolu- 
tion rule (|l|2p . therefore, defines an aggregation process 
where individuals preferably migrate, advected by v, to- 
wards the regions with higher local population. 0(s) is 
the unitary step function. A random noise £ of zero av- 
erage and correlation 

incorporates some degree of uncertainty or randomness 
in the dynamics of the system, introducing the possibility 
for the individuals to change even randomly their state. 
In a possible chemotaxis framework [23j , the system (|H2j) 
represents the motion of biological organisms that, while 
performing a random walk with diffusion coefficient D, 
secrete pheromones whose concentration turns to be pro- 
portional to the number of individuals at a given place. 
Then, the organisms following the concentration gradi- 
ent feel an attraction towards the nearest most crowded 
zones 33] . The parameter A sets the range of the veloc- 
ity excursion; the normalization factor w+(xi) + w-(xi) 
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in Eq. ([2]) makes the velocity bounded within [—A, A] im- 
plying that individuals have a limited response velocity 
as it happens in realistic cases. This maximal speed is 
reached when w+(xi) = or iy_(xj) = 0, i.e. when the 
whole population lies completely on the right/left with 
respect to Xi. The i-th individual, then, will move bal- 
listically with the absolute maximal velocity \vi\ = A to 
join the other N — 1 members independently of N. Thus 
the drift v(x{) is a nonextensive parameter defining the 
model. Equation |T]) can be cast in a dimcnsionlcss form 
with the proper rescaling of time t — > t/(a\) and space 
x — ► x/a, where only one independent parameter sur- 
vives. In the natural units (a = A = 1), the variance of 
the noise is renormalized to D — > aD/X [2(3]. Hereafter 
therefore we work with the values a = A = 1. 

In our study, we consider the effects of varying both 
the randomness and the density of individuals in order to 
see how they affects the kinetics of aggregation. For low 
noise levels (D < D c ) the system has a strong propensity 
to aggregation and almost independently of the initial 
conditions the particles form clusters that collapse each 
other to generate new larger clusters. Figure [1] provides 
some instances of the basic phenomenology of the system 
dynamics, for D = 0, D = 0.1, D = 0.5 and D = I. 

In small systems, the process continues until only one 
cluster survives with a well defined size and particle dis- 
tribution [20| | . However, in larger systems N — ► oo with 
random initial configuration, this final state is practically 
inaccessible to simulations because the transient regimes 
become extremely long. In this case the aggregation ki- 
netics is the dominating and relevant feature of the sys- 
tem behavior. 

For D > D c , the dynamics of the system loses its ten- 
dency to cluster and particles basically perform a stan- 
dard random walk inside the one-dimensional box con- 
taining them. This kind of non-equilibrium transition is 
clearly illustrated by figure [2] where we plot the evolution 
of the number of clusters N c that are present in the sys- 
tem. We define a cluster as a set of particles with mutual 
distances \xi — Xj\ below a given cutoff e that we shall 
call resolution. We see that for noise below D c — 1, N c 
fluctuates but exhibits a decreasing trend in time, signal- 
ing that clusters coalesce with each other so their average 
number decreases. 

If D > D c , after a small transient, the number of clus- 
ter fluctuates around a stationary value. A similar behav- 
ior is observed in spreading experiments where the system 
evolves from a set of spatially localized initial conditions. 
This kind of numerical simulations requires no boundary 
conditions to let the system free to occupy a broader and 
broader portion of the x-axis. The time evolution of the 
variance of the system with respect to its center of mass 



R(t) 



^2 X)<[*i(*)-*i(*)] S 



(4) 



i,j>i 



is shown in Fig. [3] for different values of D and averaged 
over independent runs. Upon rising D from zero, R(t) 
displays a crossover from a regime in which it decreases 



with time or remains bounded to a regime where spread- 
ing occurs as R(t) grows with time from the initial value 
i?(0) = <7q. This basically confirms the scenario already 
described for the number of clusters: when the dynamics 
is dominated by noise effects, the clusters become unsta- 
ble and evaporate. 

The picture emerging from the simulations on the exis- 
tence of the bifurcation point D c = 1 is affected by finite 
size effects allowing only an approximate identification of 
D c . However, the presence of a transition at the exact 
value D c = 1 can be supported by the following analy- 
sis [2(|. It is convenient to consider the Fokker-Planck 
(FP) equation associated to the Langevin equation (fT]) 
and observe that the latter can be derived formally from 
a local potential U{xi) = — ln{<p(xi)}, 



dU( Xi ) 



(5) 



where 



<f>(Xi) = ^ ex PH^ ~ Xi\) 

and D indicates the renormalized noise variance. In limit 
of large number of particles and in the mean-field like 
approach, the sum over j can be replaced by an integral 
involving the single particle density p(x,t), 



<f>{x,t) = dyp{y,t)e-\y- x \. 



(6) 



This formula clarifies that the dynamics of each walker 
can be regarded as it occurs in a time dependent mean 
potential generated by the local environment produced 
by the position of the other walkers. The integral equa- 
tion ([5]) can be transformed into a differential equation 
for <f> by a double differentiation of both members 



d 2 x 4>{x,t) - cj>{x,t) =-2p(x,t) 



(7) 



The FP equation governing the evolution of the density 
of particles obeying Eq. §5§ is 



d t p=-d x {p/<t>d x cl>)+DSP xP 



(8) 



It admits a stationary implicit solution p s {x) = 
C[(f> s (x)] 1 ^ D with C a positive integration constant de- 
termined by the normalization of p. The solution has 
to be made explicit through the substitution p s (x) into 
Eq. ©, 



-<j> s = -2C0y D 

that integrated after a multiplication for . 
bers yields 



of both mem- 



2DC 
1+D' 



(9) 



with E being another integration constant. The normal- 
izability of p s {x) and, as a consequence, of <j> s (x) for infi- 
nite systems requires that (f> vanishes enough faster when 
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FIG. 1: Trajectories of TV = 100 particles obtained by numerical integration of Eq. ([T]) with time-step h = 0.01, starting from 
uniformly distributed initial conditions in a range L = [—50,50]. The four panels, referring to parameters a = 1, A = 1 and 
noise D = 0, D = 0.1, D — 0.5 and D — 1, show the influence of the randomness on the system phenomenology. True path 
coalescence is observed only for D < 1. 
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FIG. 2: Evolutions of the number of clusters formed in a 
system with density p = 10 and different variance D of the 
noise, indicated in the inset. The threshold D c = 1 separates 
aggregation from uniform distribution of particles. 
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FIG. 3: Log-log plot of the system variance R versus time 
in spreading experiments with TV = 500 particles at different 
noise level showing the crossover from localized to delocalized 
evolutions. Each curve is the result of the average over 200 
independent runs. 



x — » ±oo implying that E — 0. Equation 0, provided 
that x is interpreted as a time coordinate, corresponds to 
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the motion of a particle with unitary mass in a potential 

v P {4>s) = -\<t> 2 s + ^ 5 <t>l +1/D ■ 

The properties of V p explain the presence or absence of 
coalescence. Indeed, when D > 1 the motion along the 
orbit E — reduces to the fixed point (f> s (x) — [p s {x) — 
0], while nonzero and localized solutions 4> s {x) can be 
found only for D < 1. This indicates the presence of a 
change in the system behaviour occurring when D goes 
through the critical point D c = 1. 

A linear stability analysis over the nonlinear FP equa- 
tion leads to the same conclusion. Let 6p{x, t) = 
exp(ikx — pt) a spatially modulated perturbation of vec- 
tor k over a uniform density background. A substitution 
into the integral gives <fi(x,t) = exp(ikx)/(l + k 2 ) 
that, in turn, plugged in the FP equation determines the 
simple condition p = k 2 (D — 1). Then the perturbation 
Sp is exponentially reabsorbed in time for D > 1 so that 
the uniform state is stable, whereas for D < 1 the pertur- 
bation is unstable signaling drastic changes in the system 
evolution. 

In the following we shall focus on the characterization 
of the aggregation dynamics for D < D c by monitoring 
and quantifying the properties of the clustering kinetics 
resulting from path coalescence. 

III. AGGREGATION DYNAMICS 

We study the dynamics of the system via simulations 
of Eqs. ([T]) and |2]) numerically integrated through stan- 
dard second order Runge Kutta algorithm for stochas- 
tic equations [jyj. All of our simulations, unless differ- 
ently indicated, started from disordered uniform configu- 
rations with particles contained in a one-dimensional box 
of a fixed size L. We used a time step h = 0.01, which 
is a fair compromise between statistical accuracy of the 
algorithm and the length of simulations. In the range 
of explored parameters we verified that a refinement of 
the time step does not affect the results. During the 
runs, when a particle close to the boundaries escapes the 
box, it is reinjected on the other side. Such reinjection 
events are rare because, in the noise range we worked, 
the drift term $2\) is dominant and the particles close 
to the boundaries move preferentially towards the bulk 
of the systems; thus border effects are really irrelevant. 
However, to check that the results are not affected by the 
re-injection through the boundaries, we repeated some of 
the simulations with fully periodic boundary conditions 
(i.e., considering system on a ring of length L). The 
different implementation of the boundary conditions did 
not alter the results. In the periodic boundary conditions 
the left and right position with respect to an individual 
located at Xi is specified by considering that each indi- 
vidual divides the ring in two halves of length L/2 each, 
accordingly, all individuals laying on the left or on the 
right of Xi can be identified without ambiguity. 



In this section, we analyze and characterize the aggre- 
gation kinetics observed in the system simulations at low 
noise regimes, where the dynamics is dominated by the 
drift term. As shown in Fig. [TJ aggregation starts by 
forming small clusters which, in turn, collide and merge 
to generate larger and larger structures. The larger clus- 
ters exhibits low coalescence rates and acts as coagulation 
centers by absorbing single particles or smaller clusters. 

A first quantitative attempt to describe the aggrega- 
tion kinetics can be made by the following argument 
that considers the very ideal case of only two clusters, 
approximated as flat particle distributions in the inter- 
vals [Xi — 6/2, Xi + 6/2] around their centers of mass 
Xi (i = 1,2). By assuming that they attract each other 
without changing too much of their structure, the aver- 
age distance between the centers of mass A = Xi — X\ 
contracts in time according to the equation (see the Ap- 
pendix A): 

A = - l + frexp ( A) < 10 > 

where b is a parameter depending on the cluster geom- 
etry, in the specific case b = 2/{l + cxp(<5/2)}; b can 
also be interpreted as a tunable parameter of a fitting 
procedure. In this ideal situation, the mean spreading, 
Eq. |gj), verifies the relation R(t) = A 2 (i)/4. In Fig. 13 
we plot the theoretical spreading R(t) obtained by the 
integration of Eq. (|10p along with the corresponding re- 
sult extracted from simulations on a system with density 
p = 10, N = 100 particles and starting from random 
configurations. We observe the formation of few clusters 
that relatively soon coalesce, A ~ 0, with a behavior 
quite in agreement with that suggested by Eq. (TIT)]) . 

However, simulations performed on larger systems to 
improve the clustering statistics revealed a rate of aggre- 
gation different from that expected from Eq. (JTUJ) that 
basically applies only to binary collisions. We, there- 
fore, investigated systematically the aggregation pro- 
cess through a set of simulations on a system involving 
N = 10 particles, where we monitored the time evolu- 
tion of the standard indicators used to quantify the de- 
gree of clustering in aggregating systems [35l [36j . Such 
indicators display a trivial behavior when the noise level 
is D > D c , as expected because uniform particle distri- 
butions there occur. 

We measured the number of clusters N c (t), their av- 
erage mass M c (t), and the distribution Pt(m) represent- 
ing the number of clusters with size m at time t. We 
also considered the average distance A(t) between neigh- 
bor clusters which provides information about the spatial 
compactness of the system. Such quantities were aver- 
aged over a set of independent runs starting from random 
particle distributions with a given density. 

In all the cases considered we observed (see Fig. [5]) the 
algebraic decay of the cluster number N c according to 

N c (t) ~ t~ z . (11) 
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FIG. 4: Time behavior of R(t) [Eq. g] for a system with 
N = 100 particles, uniformly distributed with density p — 10 
and D = 0.01. The pointed line indicates the corresponding 
quantity obtained by the numerical solution of the differential 
Eq. CEO). 
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FIG. 5: Power law decay in time of the average number clus- 
ters N c . Data refer to noise D = 0.05 and to the three densi- 
ties indicated in the inset. Lines indicate the best fit to data 
whose slopes are the exponent z. The resolution is e = 0.1 in 
all cases. 



Each cluster has a characteristic size and a typical dis- 
tribution of particles that depend on the noise and the 
other parameters in agreement with the analytical result 
for the stationary single cluster distribution derived in 
Ref. This suggests that, when sufficiently isolated, 

a cluster persists over a longer period in a sort of quasi- 
cquilibrium state. 

The results from Fig. [5J referring to systems with dif- 
ferent densities p but same noise D — 0.05, suggest an 
exponent z close to 0.2, independent of p within statisti- 
cal errors. We also observed systematically that, at each 
time t, the larger the density the lower the number of 
clusters. This means that systems with higher densities 
generate a smaller number of clusters that, due to the 
conservation of the number of particles, must contain a 
larger number of particles already in the first stages of 
their formation. The density therefore sensibly affects 
only the initial rate of cluster formation, that is, at higher 
densities clusters form faster. Later on, different densities 
reflect only in different prcfactors of the power law (fTTj) . 
This also implies that, when working at too high densi- 
ties, cluster statistics becomes worse and the power-law 
equation (fTTj) cannot be clearly detected. 

A straightforward consequence of the conservation of 
the particle number and of the decay (fTTj) is that the 
average cluster mass M c grows with the inverse power law 
M c ~ t z , as verified in our simulations (not shown for the 
sake of space). In figure [6l we plot the average distance 
between nearest neighbor clusters A directly extracted 
from the runs. Obviously, this quantity also exhibits the 
scaling 

A(i) ~ t z (12) 

with the same exponent z of the mean cluster mass. The 
presence of such a power law seems to indicate that as the 
process keeps going on and a reasonable number of larger 
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FIG. 6: Time-behavior of the average distance between clus- 
ters for different particle density, for a system with noise 
D — 0.05. The resolution is e = 0.1 in all cases, however, 
the same data processed with resolution e = 0.5 gave the 
same exponent within statistical errors. 



clusters is formed, further aggregation is suppressed and 
cluster growth becomes very slow. We also checked the 
dependence of z on the noise strength D. Fluctuations 
induced by the high noise, indeed, produce fragmentation 
or evaporation of particles from already formed clusters, 
making the influence of D on z more relevant than that 
of p. Figure [7] reports data on the variation of the slope 
[Eq. (|12p] at different D and indicates that, although 
the power law behavior remains robust with respect to 
D, until D < D Cl the z exponent lowers as D approaches 
from below the critical value D c . Therefore, in the region 
of parameter we explored, the simulation results indicate 
that the value of the exponent z weakly depends on the 
density but it is quite sensitive to the noise level D. 
To investigate in more details the aggregation process 
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FIG. 7: Power-law behavior with time of the average cluster 
distance for a system with TV = 10 4 particles and different 
noise levels D at the same density p = 1.0 and resolution 

£ = 0.1. 
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FIG. 8: Time-dependent histograms of the cluster masses col- 
lected from simulations of a system of TV = 10 4 particles, 
density p = 1.0, and D = 0.05. Lines are guides for the eye. 



we collected the histograms of the cluster masses at dif- 
ferent times during the runs. The evolution of mass 
distribution is shown in figure [8J at density p = 1 and 
D = 0.05. The shape of the distribution, at earlier times, 
rises sharply at small values of m, from zero to a max- 
imum and then decays mildly to zero at larger masses. 
The distribution are skewed and peaked around the av- 
erage and their tails are not very extended. As time 
goes on, the average increases and the distributions be- 
come broader with tails that spread over larger interval 
of masses m, indicating a clustering activity that involves 
polydispersity in cluster size. Basic theories for aggrega- 
tion phenomena (3?| suggest the scaling symmetry 



P t {m) ~ M e J{m/M c ) 



(13) 



for cluster mass distributions for processes occurring with 
an asymptotically well-defined average size M c for the 
aggregates. The constraint of the total mass (number of 



particle) conservation, 

N = ^mP ( (m), 
in 

sets the scaling exponent to the value 8 = —2 [3l|, [H, 
l38j . The rescaling of distributions in Fig. [51 according to 
Eq. (flU)) . yields the data collapse shown Fig. [5] in linear 
(top) and linear-log (bottom) plots. 

The collapse is rather satisfactory both for the bulk 
and tails of the distributions and fully consistent with 
the power-law behavior of N c (t). In fact, as an immediate 
consequence of the scaling form (fTBf , it is straightforward 
to verify that, in the continuum limit and with the change 
of variable m — > m/M c , we obtain 



N c (t) 



E 



P t (m) ~ M c (ty 



i.e., the number of clusters is proportional to the inverse 
of their average mass M c . This analysis shows that our 
model, although not directly related to general stochastic 
aggregation processes follows their typical scaling 
behavior. 



IV. CONCLUSIONS 

In this paper we considered a one-dimensional model 
that promotes aggregation of particles for a certain range 
of the parameters. The model, related to a certain class of 
problems concerning the onset of social behavior among 
individuals, is defined by a first order stochastic differen- 
tial equation. More specifically, each Brownian particle, 
representing a living organism, undergoes a drift veloc- 
ity depending nonlocally on the other particle positions. 
This drift is the response to the population unbalance 
perceived by a single individual between left and right. 
Local fluctuations occurring in the population density 
determine a displacement towards the regions with the 
highest crowding. However, a limited sensing range of 
the individuals, mimicked in the model by an exponen- 
tial weight in the distance, prevents them from reacting 
to fluctuation events occurring too far. The competition 
between far large and near, but small, aggregating centers 
makes the system dynamics nontrivial. Below a certain 
randomness threshold D c = 1, aggregation of particles in 
clusters of increasing size is observed, but at higher noise 
the organisms disperse via almost independent random 
walks suppressing social behavior. 

The aggregation process has been investigated through 
a set of simulations where we monitored the time behav- 
ior of the standard indicators used to quantify the de- 
gree of clustering in aggregating systems (35l l3o| . Such 
indicators of course display a trivial behavior when the 
noise exceeds the threshold D c because uniform parti- 
cle distributions occur. On the contrary they exhibit a 
power-law behavior at low noise regimes showing that 
this model obeys the typical scaling behavior of stochas- 
tic aggregation processes. However the relevant role of 
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aggregates, each containing half of the particles N/2. We 
assume that this situation roughly corresponds to a par- 
ticle distribution with density 

P (y, t) = p {G{S 2 /4 -(y- X x f] + Q[S 2 /4 - (y - X 2 f]} 

where po = N / {28), 0(it) indicates the unitary step func- 
tion and S represents the spatial extension of the two ag- 
gregates. Furthermore, we suppose that the intra-cluster 
dynamics is basically frozen. Therefore, each cluster 
travels as a rigid body maintaining its shape basically un- 
changed, so that the relevant time dependence lies only 
in their center mass positions X\{t) and Xz(t). This den- 
sity produces a drift velocity that can be derived from the 
integral © 



0(*i) = / dye-(y- x ^ + / 

JXt-8/2 JX 



X 2 +S/2 



dye 



-iv-Xx) 



x 2 -s/2 



then 



0(Xi) = 2(1 - e- 5/2 ) + 2sinh(<5/2)< 



m/M 



According to Eq. |5]), the velocity of cluster 1 reads 



FIG. 9: Data collapse of the time dependent distribu- 
tions shown in Fig. [8] obtained through the rescaling sug- 
gested by Eq. 1131 linear-linear (top) and linear-log (bottom) 
plots. The solid line represents the fitting curve f(x) — 
aox~ 2 exp{ai/x — C12X 2 }. 



Vi = 



sinh(<5/2) 



a Xi — X2 



0(Xi) (l-c-' 5 /2) +s i n h( ( 5/2)e x ^^ ' 
that, after some algebraic manipulations simplifies to 



the spatial fluctuations, as already noted for other pro- 
cesses Ref. [1^], does not allow one to explain the ex- 
ponents through the mean field theory defined by the 
Smoluchowsky equation for coagulation [40| . 
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APPENDIX: A 



Vi = 



1 



1 + be x ^ Xl 



(A.l) 



with b = [1 - exp(<5/2)]/sinh(5/2) = 2/[l + exp(<5/2)]. 
The velocity of the second cluster, due to the symmetry 
of the problem, has to be V2 = — Vi, so that the relative 
velocity of the two rigid clusters is — 2Vi, that is Eq. (jTUJ) 
providing that A = X2 — X\. More generally, we can 
observe that for two very localized aggregates of sizes 
Si and 82 with local density p\ and p2, respectively, the 
integrals contributing to (f> can be roughly approximated 
to pS 7 the cluster relative velocity is thus estimated as 



p 2 S 2 e 



-A 



P1S1 + p2$2e 



-A 



In this appendix we show how to derive the evolution 
equation p0[) for the centers of mass X\ and X2 of two 



The parameter b = p\&\j \p2S2) thus represents a measure 
of the relative weights of the two clusters. 
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